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ABSTRACT 

We estimate the accuracy with which various cosmological parameters can be de- 
termined from the CMB temperature and polarization data when various galactic 
unpolarized and polarized foregrounds are included and marginalized using the multi- 
frequency Wiener filtering technique. We use the specifications of the future CMB 
missions MAP and PLANCK for our study. Our results are in qualitative agreement 
with earlier results obtained without foregrounds, though the errors in most parame- 
ters are higher because of degradation of the extraction of polarization signal in the 
presence of foregrounds. 



1 INTRODUCTION 



One of the primary goals of cosmology is to accurately de- 
termine various cosmological parameters associated with the 
background FRW universe and the structure formation in 
the universe (fi, S7a, fis, ho, etc.). In recent years compelling 
theoretical arguments have emerged which suggest that the 
study of CMB anisotropics is the best hope to achieve this 
goal (Bond 1996, Knox 1995, Jungman et al. 1996, Bond 
et al. 1997, Zaldarriaga et al. 1997). On the experimen- 
tal front, two forthcoming satellite experiments MAP and 
Planck^] along with a series of ground-based and balloon- 
borne experiments on degree to sub-arcminute scales plan to 
unravel the angular power spectrum of the CMB to angular 
scales <; 1' (for details of interferometric ground-based ex- 
periments see White et al. 1997 and references therein; for 
a recent update on balloon-borne experiments see Lee et al. 
1999). It has been shown that an accurate determination 
of the CMB temperature fluctuations down to sub-degree 
scales could fix the values of nearly 10 cosmological param- 
eters with unprecedented precision (Jungman et al. 1996). 
In addition, the future satellite missions might detect the 
small, hitherto elusive signal from CMB polarization fluc- 
tuations (Bouchet et al. 1999 - hereafter Paper I - and 
references therein). The polarization data can be used to 
break degeneracy between a few parameters which are de- 
termined only in a combination using the temperature data 
alone (Zaldarriaga et al. 1997). 

One of the major difficulties in extracting the power 
spectrum of temperature and polarization fluctuations of 
the CMB is the presence of galactic and extragalactic fore- 
grounds. The extragalactic foregrounds (radio and infra-red 
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point sources, clusters, etc) will only affect small angular- 
scales (Si 10') at frequencies dominated by the CMB sig- 
nal (Toffolatti et al. 1997, 1999). The galactic foregrounds, 
on the other hand, are present at all angular scales and are 
strongest on the largest scales. They will therefore have to 
be cleaned from the future data before any definitive state- 
ments about the primary CMB signal can be made. A multi- 
frequency Wiener filtering approach was developed to study 
the implications of the presence of foregrounds for the per- 
formance of future CMB missions (Bouchet et al. 1995, 
Tegmark & Efstathiou 1996). It was shown that the pri- 
mary CMB temperature signal is much larger than the con- 
taminating foreground for all the angular scales relevant for 
future satellite missions. And therefore the performance of 
future all-sky satellite missions in extracting the CMB tem- 
perature power spectra is unlikely to be hindered by galac- 
tic foregrounds (Bouchet et al. 1995, Tegmark & Efstathiou 
1996, Gispert & Bouchet 1996, Bouchet & Gispert 1999a,b). 

In a previous paper (paper I) we extended this tech- 
nique to include polarization and temperature-polarization 
cross-correlation of foregrounds to estimate their effect on 
the extraction of CMB polarization power spectra. Our anal- 
ysis showed that the presence of foregrounds should not seri- 
ously deter the detection of _E-mode CMB polarization and 
ET cross-correlation by Planck. However, while the de- 
tection of CMB polarization will be easiest at the Doppler 
peaks of polarization fluctuations £ > 100, where it should 
help reducing the errors on the measurement of parameters 
that will already be well constrained by temperature data 
alone, the truly new information from polarization data in 
the determination of cosmological parameters is contained 
in angular scales corresponding to I < 30 (Zaldarriaga et al. 
1997). 

The polarization data helps break degeneracy between 
C2, the quadruple moment of CMB temperature fluctuation 
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and r, the optical depth to the last scattering surface (Zal- 
darriaga et al. 1997). The former gives the overall normal- 
ization of the CMB fluctuations and is fixed at the epoch 
of inflation in inflationary paradigm. The breaking of this 
degeneracy also results in a better determination of other 
inflationary parameters T/S, the ratio of tensor to scalar 
quadruple, and the tensor index nr. The optical depth to the 
last scattering surface is crucial to understanding the epoch 
of reionization in the universe. Even a value of r so small as 
0.05 leave a telltale signature in the CMB polarization fluc- 
tuations which is potentially detectable (Zaldarriaga 1997). 
The main difficulty in using polarization data to break the 
C'2-T degeneracy is that one needs to use information on 
large angular scales. The power spectra at small I is not 
only badly determined because of cosmic variance but also 
because of largely unknown level of polarized foregrounds. 

Prunet et al. (1998) attempted to model the dust po- 
larized emission from the galaxy — which is the dominant 
foreground for Planck HFI — for scales between 30' to a 
few degrees. They showed that though one might obtain 
meaningful estimates for degree scales, there can be large 
uncertainties in the large scale (£ < 50) polarized dust emis- 
sion in the galaxy. The polarized synchrotron is the other 
major galactic foreground — and it is likely to undermine the 
performance of MAP and Planck LFI. For the lack of re- 
liable data, we assumed the power spectra of polarized syn- 
chrotron to mimic that of the unpolarized component in 
Paper I. Though there remain large uncertainties on the po- 
larization foregrounds, these assumed levels of foregrounds 
combined with the Wiener filtering methods developed in 
Paper I allow us to quantify the effect of foregrounds on the 
extraction of CMB signal. In this paper, we use the methods 
developed in Paper I to ascertain the errors in the cosmo- 
logical parameter estimation. 

In the next section, we briefly review the Fisher matrix 
approach used in determining the errors on the extraction of 
cosmological parameters. We take three underlying theoret- 
ical models for our study, the rationale for which is briefly 
described in the next section. The results are presented and 
discussed in § 3. In § 4 we give our conclusions, and discuss 
the various shortcomings of our approach. 



2 FISHER MATRIX AND PARAMETER 
ESTIMATION 

Future CMB missions MAP and Planck will reach pixel 
sensitivities of ~ 30 jj,K and ~ 1.5 /J.K, respectively. This 
should allow a very precise determination of temperature 
power spectrum and a possible detection of the polariza- 
tion fluctuations (see paper I). Given the noise level and the 
underlying theoretical model, the Fisher matrix approach 
allows one to get an estimate of the errors in the estimation 
of the parameters of the underlying model. It is defined as 
an average value of the second derivatives of the logarithm 
of the Likelihood function with respect to the cosmologi- 
cal parameters, at the true parameters value (for details see 
Kendall & Stuart 1969): 



For CMB temperature and polarization data, the Fisher 
matrix can be expressed as (Tegmark et al. 1997): 



1 dC 



(2) 



where C stands for the covariance matrix of the data and 
6i correspond to cosmological parameters. The details of 
derivation of the covariance matrix and its derivatives in 
the presence of foregrounds are given in Appendix A. The 
error in the estimation of parameters is given by: 
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2.1 Underlying models 

The estimated errors on parameters will depend on the 
choice of the underlying model. We consider three models 
for our study. Though these models do not exhaust all the 
possible models and their variants, our aim is to understand 
the errors in parameter estimation for sCDM model and its 
popular variants, within the framework of generic inflation- 
ary models. We are interested in the standard parameters 
of flat FRW cosmology, h, £Ib, fii/, Oa, the reionization pa- 
rameter r, and the inflationary parameters, C2, n s , T/S, 
and tit- It is of course possible to consider a more general 
class of inflationary models which leads to a further prolifer- 
ation of inflationary parameters (Liddle 1998, Souradeep et 
al. 1998, Kinney et al. 1998, Lesgourgues et al. 1999). Wc 
also do not consider open/close universes because, as shown 
in Zaldarriaga et al. (1997), in such universes the shift in 
the angular size of the horizon at the last scattering sur- 
face leaves a very significant sign in the CMB fluctuations 
which cannot be mimicked by a change in any other pa- 
rameter, and therefore fitotai is extremely well determined 
for open/close universes. It is possible for the universe to 
be flat (or nearly flat) with contribution from both matter 
and cosmological constant, and one could attempt to mea- 
sure both these parameters from CMB data. However, the 
degeneracy between these two parameters cannot be broken 
by CMB data alone and one will have to resort to other 
measurements like observation of supernova at high-z to lift 
this degeneracy (Tegmark et al. 1998, Efstathiou & Bond 
1999). 

In addition, it is possible to include parameters like n„, 
the number of massless neutrinos, and Yf/ e , the helium frac- 
tion. However, these parameters can be better determined by 
particle physics or local observations (Jungman et al. 1996, 
Bond et al. 1997). Parameters like Q„, the contribution of 
massive neutrinos to the rest mass density in the universe, 
can be determined to a comparable accuracy by the data of 
future galaxy surveys like SDSS (Hu et al. 1998). 

In this paper, we take only CMB data for our study 
and do not include priors from other measurements like fu- 
ture Galaxy surveys or high-z Supernova results. The three 
models we consider are: 

1. sCDM model with r = 0.1. The rather large value of r 
is taken to bring out the effects of polarization data. 

2. Tilted CDM model with r = 0.1, n s = 0.9, n T = -0.1, 
and T/S = 0.7. Note that T/S = -7n T , which is one of the 
predictions of slow-roll inflation (Starobinsky 1985, Liddle 
& Lyth 1992). 
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3. Model 2 with f2„ = 0.3 with two light, massive neutri- 
nos. 



3 RESULTS 

We use the results of Paper I (the values of various terms 
in the covariance matrix as defined in Appendix A) for the 
specification of the future satellite missions. Our results are 
shown in Tables [j], ^| ^ for the three underlying models. 

It should be noted that we fix the value of C2 = 
Ci + Cf = 796 {fiKf for all the models. Only for the 
sCDM model does it correspond to COBE normalization. In 
the models with tensor contribution, the COBE normalized 
CMB signal is larger than the signal for our normalization 
by a factor of ~ f.5. 

To assess the reliability of our code we computed the 
errors in the 6-parameters sCDM model of Zaldarriaga et 
al. 1997 with their instrumental specifications and without 
foregrounds, and compared our results with both theirs and 
those obtained by Eisenstein et al. 1999. Our results are 
comparable to those of Eisenstein el al. 1999 for most pa- 
rameters, with the exception of r where the error we find is 
noticeably bigger. We think that this discrepency is related 
to the way we compute the derivatives of the spectra with 
respect to the parameters (see appendix A). 

The results for sCDM model are shown in Table [j]. For 
comparison, results for the best channel of each experiment 
without foregrounds are also shown. As is clearly seen, the 
performance of Wiener filtering matches the best channel 
case for all the experiments. In Paper I we showed that the 
Wiener filtering performance in extracting the temperature 
power spectra lies between the expected performances of the 
best channel and the combined sensitivity of all channel for 
each experiment, at least for the specific foreground models 
considered. As the temperature data alone gives a fair idea 
of the errors on most of the parameters our results could be 
anticipated from conclusions of Paper I. However, the errors 
of C2 and r are mostly determined by the polarization data. 
In paper I we showed that the extraction of polarization 
power spectra is degraded as compared to the cases with no 
foreground. Our results in this paper suggest that it should 
not be too much of a deterrent in determining cosmologi- 
cal parameters. It is also important to note that the present 
results for the best channel case are comparable to Wiener 
filtering case. This means that i) the other channels can ef- 
fectively be used to clean the best channel ii) the presence of 
foregrounds do not introduce additional degeneracies which 
are absent when the data is assumed to contain only CMB 
and pixel noise. 

In Table ^, the expected errors are shown for a model 
which includes contribution from tensor modes. One of the 
aims of studying this model is to establish how well the infla- 
tionary parameters can be determined. In comparison with 
the sCDM case, the errors on all the standard parameters 
are bigger. This is because the additional parameter T/S 
allows one to fix the normalization more freely, thereby in- 
troducing additional degeneracies (Zaldarriaga et al. 1997). 
The errors on parameters like C2, h and Q,b are higher than 
for similar models considered by Zaldarriaga et al. (1997). 
This is partly due to our choice of normalization which gives 
smaller signal. However, it also reflects the degradation of 



the polarization power spectra extraction in the presence of 
foregrounds. Other parameters like T/S, t, and nr are bet- 
ter determined than the results of Zaldarriaga et al. (1997), 
but it is mostly owing to the fact that we take larger in- 
put values for r and T/S. We also show the effect of in- 
cluding very small signal from B-mode polarization. As is 
seen, it results in a better determination of most parame- 
ters, especially the inflationary parameters. Though the B- 
mode signal is much smaller compared to _E-mode signal, 
and is generally below the pixel noise except for a small 
range of modes for i Sa 100, its very presence indicates ten- 
sor modes in inflationary paradigm. Also, the degradation of 
the extraction of this signal in the presence of foreground is 
smaller than for the iJ-mode signal (Paper I). Therefore, it 
can make a difference in the estimation of parameters. Our 
results show that the consistency condition of slow-roll in- 
flation, T/S — ~7n,T, can be checked by future missions (it 
should be noted here that this relation was imposed only in 
the fiducial model, but excursions of both parameters were 
considered independently). Planck can extract both these 
parameters with la errors S3 50%. However, it should be 
kept in mind that our results are more optimistic than the 
results of Zaldarriaga et al. (1997) because of our choice of 
input model. 

The results of adding another parameter Q. v in the 
model above are shown in Table |S| f2„ can be determined 
to an accuracy ;S 10% with Planck, though it will be very 
difficult for MAP to determine it. Note that errors on other 
parameters have not changed much by the addition of this 
parameter, which suggests that no new degeneracies have 
cropped up. However, degeneracies between various param- 
eters depend very sensitively on the choice of input model. 
For instance, if the new parameter Q. v was added with the 
input value f2„ = 0, it would have substantially worsened 
the estimation of almost all parameters, especially the in- 
flationary parameters. For all the three models considered 
here, we took = 0. A finite value of Qa results in a better 
estimation of all the parameters of FRW model as well as 
Qa (Zaldarriaga et al. 1997). 



4 CONCLUSIONS AND SUMMARY 

In this paper we estimated the effect of foregrounds on the 
determination of cosmological parameters. The most im- 
portant result is that although the presence of foregrounds 
somewhat worsens the parameter estimation by degrading 
the detection of polarization signal, it does not give rise to 
any severe degeneracies not already present in the CMB data 
(CMB signal and pixel noise) without foregrounds. It needs 
to be further confirmed with a detailed study of the Likeli- 
hood function in the multiple parameter space. 

Any analysis such as ours can only give a qualitative 
idea on the accuracy of parameter estimation. This is largely 
because of its strong dependence on the input model (Zal- 
darriaga et al. 1997). In addition, there are great uncer- 
tainties in the assumed level of foregrounds which we take 
in the Wiener filtering analysis of paper I. Moreover, the 
foreground characteristics (power spectra, frequency depen- 
dence) should be determined from the data as well, and this 
adds uncertainty to the determination of the cosmological 
parameters (see Knox 1999). It should be noticed at this 
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point that, if Wiener filtering assumes some frequency de- 
pendences as well as power spectra for the CMB and the 
foregrounds, it also gives a measure of the error on the esti- 
mation of the power spectra of foregrounds from the filtered 
data, see Paper I. ' 

Still our results suggest that the primary obstacles for 
high precision CMB measurements will rather stem from 
systematic errors and inaccuracies in calibration, baseline 
drifts, determinations of far side lobes or estimates of filter 
transmissions. . . all of which are of course not included in 
this analysis. Furthermore, a Fisher matrix analysis leads to 
the smallest possible error bars, which can only be degraded 
by inaccuracies in devising the Wiener filters (by using ap- 
proximate power spectra and spectral behaviours). 

The next step will be to directly analyse simulated 
mega-pixel multi-frequency CMB maps relevant to future 
experiments. However, such an analysis is an extremely dif- 
ficult (if not intractable) numerical problem (for recent at- 
tempts see Muciaccia et al. 1997, Oh et al. 1999, Borrill 
1999). In light of this, our results should be regarded as 
a first attempt on the problem of parameter estimation in 
the presence of foregrounds, which give a qualitative idea 
of the expected accuracy in parameter estimation till the 
analyses of multi-frequency CMB datasets become possible. 
Since the submission of this work, more detailed analysis of 
the effect of foregrounds have been investigated by Tegmark 
et al. 1999. Their results are similar to ours, with maybe 
slightly higher foreground residuals as they allowed some 
scatter in the frequency dependence of foregrounds. 



5 APPENDIX A 

In this section we briefly recapitulate the discussion of paper 
I, and derive the covariance matrix of CMB data and its 
derivatives. The observed CMB data at multiple frequencies 
can be expressed in multipole space as: 



„' ( x p x i>) 



(6) 



yl(l,m) = Al J p (l,m)x p (l,m) +6*(Z,m) 



(4) 



where x p is the underlying signal for process p and "field" 
j (i.e. temperature or (E,B) polarization modes), and v is 
a frequency channel index. In the Wiener filtering method, 
one considers a linear relation between the true, underlying 
signal, x p and the linearly optimal estimator of the signal, 
Xp . 

4 = w;iyi. (5) 

Eqs. (Q) and (|s|) can be used to write the estimated power 
spectrum as: 

(4x pl ) = (WA)™„(WA)»,„(^^„,) 

+w; l v w^,{b l v bz,) 

t It sould also be noticed that any spatial change of the frequency 
indexes (for dust or synchrotron) should correspond to special 
astrophysical regions (molecular clouds, supernovae remnants. . . ) 
and that an analysis where any such spatial change of index is 
simply incorporated as an additional "noise" term would lead 
to pessimistic results. This would rather point out that a global 
analysis (in "Fourier" modes) is insufficient to take this effect 
properly into account. 



where the last equality comes from the equation defining the 
Wiener filter (see Eq. 6 of Paper I). The covariance of the 
filtered data can then be written as: 



(QVCti Q\ 2 Ctei 
Q 12 Ctei Q 22 Cei 







Qi^Cbi 



(7) 



For computing the Fisher matrix we also need to com- 
pute the derivative of the covariance with respect to cosmo- 
logical parameters: 

dC £ _ _dC e dC e (X) 



X=T,E,ET,B 



dC e (X) ddi 



(8) 



The derivative of the covariance matrix with respect to var- 
ious power spectra can be written using Eq. ([]). These 
derivatives, it should be borne in mind, are with respect 
to the input power spectra used in estimating the Fisher 
matrix and not the power spectra used in constructing the 
Wiener filters, which, therefore, are invariant under this 
change. These derivatives can be readily calculated: 



d ( x p x p) = (W^A 11 ) 2 



2 x (w£a%w£a%) 



dCj E 

9( > x p x p) = (w 12 A 22 ) 2 

Ql-E E\ 

°\ x p x p> = ( W 22 A 22 ) 2 

yE \ yv pv ^Vp ) 



p 

d{x^xf) y (w 22 A 22 W 21 A 11 i 

Q(jT E \ pu -^vp ' * pu -^-upj 

di^Xp Xp ^ 

acj 

3(Xp Xp } 

d{Xp Xp ) 
~dCf~ 



= (W^All) 2 

= (W^A^W^A 11 ) 

\ * ' pu J L isp ' * pv J *-vp } 



(9) 
(10) 

(11) 
(12) 
(13) 
(14) 
(15) 



(W 11 A U W 22 'A 22 ) + (W 12 A 22 W 21 A X1 V16) 

\* ' pu J *-vp f * pv J *-vp ) 1 V V v vp * * pv J vp) \ ^ I 



\^pu A-vp ^^pv ^-up) 



(17) 
(18) 



Theoretical power spectra are calculated using the CMB 
Boltzmann code CMBFAST (Seljak & Zaldarriaga 1996). 
Derivatives with respect to cosmological parameters are cal- 
culated numerically using a variant of df ridr routine of nu- 
merical recipes (Press et al. 1992). We notice that a 5% 
step in most parameters gives stable results. The only ex- 
ception is derivative of i?-mode power spectra with respect 
to t when r < 0.05 for £ < 20. This numerical instability is 
expected as a small change in this parameter when the input 
value of r is very small can cause appreciable change in the 
.E-mode power spectra at small I. However, the numerical 
differentiation is quite stable for r <; 0.05. 
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Table 1. Errors on parameters for sCDM model with r = 0.1. 
Also shown are the corresponding errors for the best channel of 
each experiment. 



Parameters 


C<2 


h 
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n s 


Model 


/yo^ii ) 


U.O 


U.Uo 


U.U 


U. 1 


l.U 


wiener lainckj 


A oz, 
/o 


l.OO /o 


z.o /o 


u.uoy 


4.D /0 


U.o4 /o 


Best Channel (Planck) 


2.1 % 


1.06 % 


1.82 % 


0.03 


3.74% 


0.3 % 


Wiener (HFI) 


2.48 % 


1.39 % 


2.37 % 


0.04 


5.75 % 


0.35 % 


Best channel (HFI) 


2.1 % 


1.06 % 


1.83 % 


0.03 


3.75 % 


0.3 % 


Wiener (LFI) 


3.81 % 


2.26 % 


3.72 % 


0.067 


10.3 % 


0.54 % 


Best Channel (LFI) 


3.6 % 


2 % 


3.3 % 


0.057 


6.5 % 


0.51 % 


Wiener (MAP) 


4.9 % 


4 % 


8.9 % 


0.12 


45.6 % 


1.65 % 


Best Channel(MAP) 


6.3 % 


4.5 % 


10.7 % 


0.13 


43.5 % 


1.76 % 



Table 2. Errors on parameters for a model with tensor contribu- 
tion with or without the inclusion of B-mode polarization 



Parameters 


c 2 


h 


n 6 


n A 


r 


«S 




T/S 




Model 


796(pK) 2 


0.5 


0.05 


0.0 


0.1 


0.9 


-0.1 


0.7 




Wiener (Planck) 


8.7 % 


1.6 % 


2.7 % 


0.045 


5.5% 


0.46 % 


81 % 


22.4 % 




+ B-modes (Planck) 


6.5 % 


1.55 % 


2.64 % 


0.044 


4.8% 


0.43 % 


57.1 % 


17.5 % 




Wiener (HFI) 


9.4 % 


1.63 % 


2.8 % 


0.05 


7 % 


0.47 % 


87 % 


24.1 % 




+ B-modcs (HFI) 


7.7 % 


1.6 % 


2.7 % 


0.05 


6 % 


0.45 % 


70 % 


20.6 % 




Wiener (LFI) 


9.8 % 


5.3 % 


8.6 % 


0.15 


11.3 % 


1.65 % 


91.6 % 


32.4 % 




+ B-modcs(LFI) 


9 % 


4.6 % 


7.5 % 


0.13 


9.6 % 


1.42 % 


83 % 


28.2 % 




Wiener (MAP) 


12.3 % 


22.3 % 


40 % 


0.67 


52.5 % 


7.5 % 


91 % 


91 % 




+ B-modcs (MAP) 


12 % 


20 % 


36.5 % 


0.60 


46 % 


7 % 


90.4 % 


81.6 % 




Table 3. Parameter estimation with Q 


















Parameters 


C 2 


/i 


n b 


n A 


T 


«s 


nj" 


T/S 




Model 


796(^) 2 


0.5 


0.05 


0.0 


0.1 


0.9 


-0.1 


0.7 


0.3 


Wiener (Planck) 


8.8 % 


1.6 % 


3 % 


0.045 


4.6% 


0.95 % 


82 % 


24 % 


10.5 % 


+ B-modcs (Planck) 


6.4 % 


1.45 % 


2.65 % 


0.04 


4.3% 


0.81 % 


55 % 


18 % 


9.45 % 


Wiener (HFI) 


9.6 % 


1.66 % 


3.1 % 


0.045 


6 % 


0.97 % 


89 % 


26 % 


10.8 % 


+ B-modcs (HFI) 


7.75 % 


1.55 % 


2.9 % 


0.42 


5.5 % 


0.9 % 


71 % 


22.6 % 


10 % 


Wiener (LFI) 


9.2 % 


3.6 % 




0.098 


11 % 


2 % 


84 % 


31 % 


26.5 % 


+ B-modcs(LFI) 


9.2 % 


3.6 % 


6 % 


0.097 


10.8 % 


2 % 


84.8 % 


31 % 


26.5 % 


Wiener (MAP) 


12.3 % 


18.4 % 


27 % 


0.69 


67.1 % 


7.5 % 


115 % 


75 % 


201 % 


+ B-modcs (MAP) 


12 % 


17 % 


25.4 % 


0.65 


62 % 


7.1 % 


115 % 


70 % 


200 % 
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